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Abstract 

The collective interactions of many partons in the first stage of 
the collisions is the usual accepted explanation of the sizable elliptical 
flow. The clustering of color sources provides a framework of par- 
tonic interactions. In this scheme, we show a reasonable agreement 
with RHIC data for px < 1.5 GeV/c in both the dependence of V2 
on transverse momentum and in the shape of the nuclear modified 
factor on the azimuthal angle for different centralities. We show the 
predictions at LHC energies for Pb-Pb. In the case of proton-proton 
collisions a sizable V2 is obtained at this energy. 



1 Introduction 

A major breakthrough was the discovery by RHIC experiments of a large 
elliptic flow t> 2 [lj|2j|3j|4j [5j[6j[7j[8j[9j. A non vanishing anisotropic flow 
exist only if the particles measured in the final state depend not only on the 
physical conditions realized locally at their production point, but if particle 
production does also depend on the global event geometry. In a relativistic 
local theory, this non local information can only emerge as a collective effect, 
requiring interactions between the relevant degrees of freedom, localized at 
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different points of the collision region. In this sense, anisotropic flow is a 
particularly unambiguous and strong manifestation of collective dynamics 
in heavy ion collisions [TUJ. The large elliptic flow v-i can be qualitatively 
understood as follows: In a collision at high energy, the spectators are fast 
enough to free the way, leaving behind at mid- rapidity an almond shaped 
azimuthally asymmetric region of QCD matter. This spatial asymmetry 
implies unequal pressure gradients in the transverse plane, with a larger 
density gradient in the reaction plane (in-plane) perpendicular to it. As 
a consequence of the subsequent multiple interactions between degrees of 
freedom, this spatial asymmetry leads to an anisotropy in momentum space; 
the final particle transverse momentum are more likely to be in-plane than 
out-plane, hence v 2 > as it was predicted [TT] . 

The experimental results show that the elliptical flow first increases with 
Pt and then levels off around px ~ 2 GeV. Also, the p?- dependent v% of iden- 
tified hadrons shows mass ordering at small px and displays a constituent- 
quark counting rule. On the other hand v<i normalized to the eccentricity 
scales on the density of produced charged particles to the overlapping area 
of the collision. To understand these experimental results, many theoreti- 
cal approaches have been used. These include hydrodynamics models [T2] . 
transport models [13], and hybrid models which combines both hydrody- 
namic and transport models [H]. Hydrodynamics models usually gives the 
largest elliptical flow. 

In transport models, including the parton cascade [T5], it is obtained 
a sizable v-i only if it is used a very large parton cross section. On the 
other hand, transport or hybrid models based on string degrees of freedom 
in general gives a smaller elliptic flow than the observed at RHIC |16j . 

The inclusion of interaction among strings via string fusion as in the 
AMPT model[T7|, in dual parton model [IB] in the quark gluon string model 
[T9] or in Reggeon theory [20j yield a large elliptic flow. In this paper, we 
would like to study this point further, using percolation of strings. In this 
model, is introduced the interaction between many degrees of freedom. In 
fact the overlapping of the strings forming color clusters, can be seen as an 
interaction among the partons of the projectile and target, from which the 
strings are stretched. The recombination of the flavour and color of these 
partons give rise to the flavours and color content of the resulting cluster. Due 
to that, the fragmentation of the clusters is different of the fragmentation of 
the initial strings, corresponding to a higher color and also different flavour. 
This fact, gives rise to a suppression of the multiplicities, a moderate increase 
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of the transverse momentum and an increase of the ratio baryon/meson with 
centrality. Also it provide us with an state which has strongly interacted in 
the first initial times. The percolation of strings incorporates in a dynamical 
way two of the main features of AA collisions, namely recombination and 
the formation of a strong interacting partonic system. We will compare 
the results obtained for the elliptic flow with experimental RHIC data on 
AA collisions, including the dependence of the nuclear modified factor on 
the azimuthal angles for different transverse momentum. We will show our 
predictions for LHC energies, which gives a slightly larger elliptic flow for 
all particles. A proton can be considered as an extended object, although of 
smaller size than a nucleus. In proton- proton collisions at LHC energies a 
high number of strings are stretched between the partons of similar effects 
to those seen at RHIC experiments could be expected. 

At first sight the difference between px spectra observed in pp and AA 
collisions at RHIC energy would suggest little flow, a conclusion in line with 
naive expectations. However, it has been shown [2TJ that the effects of energy 
and momentum conservation actually dominate the observed systematic and 
that pp collisions may be much more similar to heavy ion collisions than 
generally thought. In the framework of percolation of strings we derived an 
approximate analytical formula for the pt distributions, valid for AA and pp 
collisions. Using this formula we compute the elliptic flow for pp collisions, 
at LHC energy, obtaining a sizable size. This is not unexpected because the 
string density at LHC energy in pp collisions is close to the critical percolation 
density and therefore there are many strings which overlap each other, giving 
rise to high density partonic medium. 

The plan of the paper is as follows: In the next section, we introduce the 
percolation of strings, computing the elliptic flow in section 3. In section 4 we 
discuss our results comparing with the RHIC data and show our predictions 
at LHC, for AA and pp collisions. Finally we present our conclusions. 

2 The string percolation model 

In the string percolation model [22] [23] [21], multiparticle production is 
described in terms of color strings stretched between the partons of the pro- 
jectile and the target. These strings decay into new ones by qq or qq + qq 
pair production and subsequently hadronize to produce hadrons. Due to the 
confinement, the color of these strings is confined to small area in transverse 
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space Si = ttTq with r ±2 0.2 — 0.3 fm. These values comes from lattice results 
on bilocal correlators [25] and from the method of field correlators [26] [27] 
and corresponds to the correlation length of the QCD vacuum. With in- 
creasing energy and/or atomic number of the colliding particles, the number 
of strings, N s , grows and they start to overlap forming clusters, very much 
like disks in two dimensional percolation theory. At a certain critical density, 
a macroscopical cluster appears, which marks the percolation phase transi- 
tion. This density corresponds to the value p c = 1.18 — 1.5, (depending on 
the type of employed profile functions of the nucleus, homogeneous or Wood- 
Saxon [28]) where p = N s Si/Sa and Sa stands for the overlapping area of 
the colliding objects. The percolation theory governs the geometrical pattern 
of the string clustering. Its observable implications, however, require the in- 
troduction of some dynamics in order to describe the behavior of the cluster 
formed by several overlapping strings. We assume that a cluster of n strings 
behaves as a single string with an energy-momentum that corresponds to the 
sum of the energy-momenta of the individual strings and with a higher color 
field, corresponding to the vectorial sum of the color field of each individual 
string. In this way [2H], we can compute the mean multiplicity p n and the 
mean transverse momentum squared < p\ > n of the particles produced by a 
cluster. 



nS n 2 nSi o , . 

< Hn >= y -g- < Hi > and < Ptu >= y — < Pn >, I 1 ) 

where < pi > and < "Pri > are t ne mean multiplicity and mean p\ of particles 
produced by a single string. We take Si constant equal to a disc of radius 
r . S n corresponds to the total area occupied by n discs, which of course 
can be different for different configurations even if the clusters have the same 
number of strings. Notice that if the strings are just touching each other 
S n = nSi and the strings act independently to each other. On the contrary 
if they fully overlap S n = Si it is obtained the largest suppression of the 
multiplicities and the largest increase of the transverse momentum. In the 
limit of high density, one obtains 



Si P 1 / r-\\ 

< >= — = fW (2) 

and the equations (1) transform into the analytical ones. 

<fi>= N s F(p) < 111 >, and < p 2 T >= J\ (3) 
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with F(p) = y 1 - p P . In principle, the equality (3) has been obtained in 
the high energy limit, what means close or above, the critical density, where 
the strings covers a large fraction of the total available area. However the 
equations (3) give also the right behavior in the limit of low density, in fact 
if p tends to zero, F(p) is one and therefore there is not any change from the 
result of the superposition of individual strings. In the mid rapidity region, 

4 

the number of strings is proportional to the number of collisions, N co u ~ iV| 

2 

being Na the number of participants. Therefore p ~ iV| and p ~ In 
other words, the multiplicity per participant becomes independent of Na, i.e. 

4 

saturates. Outside the mid rapidity, iV s is proportional to Na instead of iV| . 

i 

Therefore, there is an additional suppression factor iVf compared to central 
rapidity. 

If we are interested in a specific kind of particle i, we will use equations 
(1) and (3) pu, < p^i < A*n >% and < p\ n >i for the corresponding 
quantities. 

In order to compute the multiplicities we must know, N s and p± (for a 
fixed centrality, knowing N s we deduce the density p ). Up to RHIC ener- 
gies, N s in the central rapidity region is approximately twice the number of 
collisions, N co u. However N s is larger than N co u at RHIC and LHC energies. 

Using equation (3) the multiplicities are computed using a Monte-Carlo 
code based on the Quark- Gluon String model to generate the strings |30j. 
Each string is produced at an identified impact parameter. From this, know- 
ing the transverse area of each string, we identify all the clusters formed in 
each collision and subsequently compute for each of them its multiplicity in 
units of p\ the value of p\ was fixed by the comparison of our results with 
SPS data for Pb-Pb central collisions. Our results are in agrement with SPS 
and RHIC multiplicity data [3Tj . Using the first of equations(l) we obtain 
very similar results [32]. Must of the reasonable string models [30] [33] [34] 
[35J[36J obtained similar results for N s . This fact gives us confidence is our 
values. Notice that sometimes, even in experimental analysis, N co u is ob- 
tained from the optical approximation of the Glauber model without taking 
into account the energy-momentum conservation. This conservation reduces 
N co u. In the evaluation of the elliptical flow, we need the values of dN/dy at 
central rapidity for the different collisions and centralities, and values of p. 
We take both from our previous results [31 J [32J . 

Concerning the transverse momentum distribution, one needs the distri- 
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bution g(x,px) for each cluster, and the mean squared transverse momen- 
tum distribution of the clusters W(x), where x is the inverse of the mean 
of the transverse momentum of each cluster which is related to the clus- 
ter size by equation (3). For g(x,pr) we assume the Schwinger formula, 
9{ x iPt) = exp(—PT x ) [23] [37] used also as a first approximation for the 
fragmentation of Lund string [22] • For the weight function, we assume 
the gamma distribution W(x) = 7 ^L — e~ fca: ^.[55] [39] with 7 = and 



7(fc) it?£*i ,viu " / " <x> 

is proportional to the width of the distribution, depending 



(<x 2 >-<a:> 2 ) 



<x> z 



on the density of strings r\. Therefore the transverse momentum distribution 
f(pT,y) of the particle i is 

f(PT,y) = -QTfi = f™dxW(x)g(p T ,x) 

= 1 F{p) ! (4 ) 

dy k <p 2 T >u {P) (l + F(p)p 2 T /k<p 2 T > li f {) 

The equation (4) is valid for all types of collisions and high densities. The 
equations (4) can been seen as superposition of color sources, clusters, where 
i fixes the transverse momentum fluctuations. At small density of strings, 
there will not be overlapping of strings, the strings are isolated and k — > 00. 
When p increases, there are some overlapping of strings, forming clusters 
and therefore k decreases. The minimum of k will be reached when the 
fluctuations reaches its maximum. Above this point, increasing p, these 
fluctuations decrease and k increases. In the limit of p — > 00, there is only one 
cluster formed by all the produced strings. Again, there are no fluctuations 
and k tends to infinity. 

The quantitative dependence of k on p was obtained in [SB] from the 
comparison of equation (4) with RHIC AA data at different centralities. We 
use here such a dependence. The peripheral Au-Au collisions at y/s = 200 
GeV, corresponds to the p values slightly above the minimum of k [38J . In 
the range of the used p values, k is a smooth function on p. We will show 
that for pt < 1.5 GeV/c our t> 2 does not depend on k. The formula (4) 
is able to describe the transverse momentum distribution for mesons in AA 
collisions at all centralities and energies up to moderate px, Pt < 4—5 GeV/ c. 
In particular they describe rightly the Cronin effect and the suppression of 
moderate pt particles, as it was shown in reference [38] . 

We do not claim to describe the data at high p?- It is well known that 
jet quenching is the mechanism responsible for the high pt suppression. This 
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phenomena is not included in our formula (4) which was obtained assuming 
a single exponential for the decay of a cluster without a power like tail. Our 
formula can be seen as a way to join smoothly the low and moderate pr 
region with the high px region. Indeed, the high pt suppression implies by 
continuity a suppression of the highest pt values of the intermediate region 
which are described by formula (3). 

The universal formula (4) must be considered as an analytical approxi- 
mation to a process that consists on the fragmentation of clusters of strings 
and their eventual decay via successive cluster breaking via Schwinger mech- 
anism [3D]. The formula (4) is not able to explain the differences between 
antibaryons (baryons) and mesons. In fact, the only parameter that is dif- 
ferent for them in (4) is the mean of the squared of transverse momentum 
of pions and protons produced by a single string < p\ >\ n and < p\ >\ p , 
respectively. This only causes a shift on the maximum of Raa but keeps the 
same height at the maximum contrary to what one observed. However, in 
the fragmentation of a cluster formed by several strings, the enhancement of 
production of antibaryons (baryons) over mesons is not only due to a mass 
effect corresponding to a higher tension due to a higher density of the cluster 
(the factor F(p) in front of p\ in (4)). In fact, the color and flavour proper- 
ties of a cluster follow from the corresponding properties of their individual 
strings. A cluster composed by several quark-antiquarks (q — q) strings be- 
haves like an string with a color Q and flavour composed of the flavour of 
the individual strings, as a result, we obtain clusters with higher color and 
different flavour ends. For the fragmentation of a cluster we consider the 
creation of a pair of parton complex as QQ [3DJ , after the decay, the two new 
QQ strings are treated in the same manner and therefore decay into more 
QQ strings, until they come to objects with masses comparable to hadron 
masses, which are identified with observable hadrons by combining with them 
the produced flavour with statistical weights. In this way, the production of 
antibaryons (baryons) is enhanced with the number of strings of the cluster. 
The additional antiquarks (quarks) required to form an antibaryon (baryon) 
are provided by the antiquarks (quarks) of the overlapping strings that form 
the cluster [40] [41]. Therefore recombination or coalescence mechanism [32] 
[43] are incorporated in a dynamical way. In order to take this into ac- 
count keeping our close formula we will do the following approximation. For 
antibaryons (baryons) we consider instead of the equation (3), 

Hb = Nl +a F( PB )fi 1B (5) 
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fitting the parameter a to reproduce the experimental dependence of the px 
integrated p spectra with centrality [3H]- The obtained value is a = 0.09. this 
means that the density p must be replaced by p§ = Nfp The antibaryons 
(baryons) probe a higher density than mesons for the same energy and type 
of collision. The equation (4) with the described modification for baryons 
provide us the input for the evaluation of the elliptic flow. 



3 Elliptic Flow. 

In saturation initial state models as the color Glass Condensate [H] [45] or 
the Percolation String Model [22] [23] color sources radially emit particles. 
The azimuthal asymmetry has, in this case, the origin in the deviation from 
the circle of the impact parameter projected interacting overlap region in 
non central, b > 0, collisions. The situation is schematically represented in 
Fig. 1. If one measures production within fixed azimuthal angles, one has to 
know the source density as a function of R^. If no azimuthal determinations 
are made one can adopt a circle of radius R such that 



ttR 2 _ 1 
~T~ ~ 2 



R'dip (6) 



or 

R 2 =< K > (7) 

Note that, as in the experiment, we are using for the range of variation of 
(p, < (p < |. As the number of color sources is the same in the interior of 
the circle and in the ellipsoid and as the directional density is related to the 
average density by the formula 

P V = P(#) 2 (8) 

one expects higher < p\ > and larger particle numbers for production in the 
reaction plane, (p — 0. In order to make this point clear let us consider a 
model, simpler than (4), which is the Schwinger model with percolation [37J: 

2 dN -F(p)^— 

-^~e <^ 9 

7T dp j, 
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for the azimuthal integrated cross section and 
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, 9 , ~ e <p t>i (10) 
dp^dip 

for a direction <p. From (9) and (10) one sees that the ratio of the px inte- 
grated quantities is 

(10) < p\ > y F(p) 
(9) <^> F(p,) 

which for small <p, as -F(p) is decreasing function of p, is larger than 1. 
The same result works for Eq. (4). The relevant directional factor is the 
factor F(p)p T or F(p (/ ,)p|,. Let us define, for a given centrality, number of 
participants or p, 

4 Us/(rw ' l4) (12) 

and 

^ s/ (^) (13) 

We will treat the azimuthal distribution as perturbation around the average, 
isotropic distribution. Such expansion, as p ( is proportional) to and p^ 
to -jj-, can be written as: 



din f(p 2 T ,R 2 ) ( 2 2 . 
! + ^ (^-^) 



(14) 



Note that (4), because of (6), satisfies the normalization condition, 

[ W/2 f(F(p^),p 2 T ) d<p = f {F(p),pl) (15) 

Additional terms in the expansion will, in general, violate the normalization 
condition (15). As a decrease in R 2 corresponds to an increase in source 
density and an increase in particle production one expects, in (16) 

< 

dn z 

with the consequence, again see (16), that production is maximal for ip ~ 0. 
Note that the second term in the right hand side of (16), changes sign at 
= R 2 . 
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Figure 1: 



We shall next write expressions for relevant measurable quantities as mod- 
ification factors, Raa and v-i- With the usual definition for R AA (p T ), we 
define, 



AipR AA (p%) 



and 



V2 (Pt) 



2 



7T 



dip cos(2ip) 



\ + d In f AA (F(p),p 2 T ) , r2 _ r2 , 
dR 2 ^ v ' 



(17) 



;is) 



Note that (16), using (4), becomes 



f(p 2 T ,Rl)=f(p 2 T,R 2 ) 



I , „ „, ,ox Pt 1 



i? 2 F(p) 



2p?l + 



(19) 

From the equations (19), (20), (21) we can compute our main observables to 
compare with the experimental data. The approximation done is valid when 
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the second term in (21) is small, what means not hight pt- Notice that due to 
the location in impact parameter space of the strings that form clusters there 
is an initial spatial azimuthal asymmetry which is translated into transverse 
momentum due to interactions of strings included in the formula for the pr 
distribution via the factor F(p). The strength of the interaction depends 
locally on the string density and therefore is the origin of the azimuthal 
asymmetry. 

4 Results and Discussion 

In order to compute the elliptic flow, we need the values of N s , p, k(p) and 
< Pt > 1vd < Pt >ifc> < P% >ip- As we mentioned before, for each type 
of collision, energy and centrality, we obtain the number of strings using 
a Monte-Carlo code based on the Quark Gluon string model [30], hence 
we compute p and k(p). The dependence of k on p is very smooth, our 
results for p T < 1.0 GeV/c are independent of k according to formula (18) 
and (19). In table 1, we show the value of N s , p, k(p) for different types of 
collisions, centralities and energies. We include also the values of the charged 
multiplicity at central rapidity ^ which we will use below. They are taken 
from previous work done in the framework of percolation [31] [32] |37j[46j. 
We use the values < p\ > l7T = 0.0Q(GeV/c) 2 , < p\ > lk = O.U(GeV/c) 2 , < 
Pt >ip= 0.30(GeV/c) 2 , obtained previously from the comparison of formula 
(4) with the experimental data of the dependence of the mean transverse 
momentum of the different particles on the multiplicity [37] [38J. 

Once these parameters are fixed then the elliptic flow for a certain type of 
collision, centrality and energy is known because it only depends on p. Notice, 
that we could consider < p\ > 7r , < p\ > lk , < p\ > 1)5 as free parameters doing 
a fit to the experimental data. However, our aim is not to obtain a perfect 
agreement with data. We are aware of some simplifications and limitation 
of our framework, therefore we just check whether the general trend of the 
data is reproduced in our model. 

In fig. 2 we show our results for minimum bias Au-Au collisions for 
pions, kaons and antiprotons compared with the experimental data |47| . A 
reasonable agreement is obtained for pt < 1 GeV/c, reproducing the mass 
ordering. In fact, at low pt, the elliptic flow is proportional to p T and the 
proportionality coefficient is a . As < p T1 >tt<< p T i >k<< Pti >p 
t>2 is larger for pions than for kaons and antiprotons. Our evaluation was 
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Figure 2: Elliptic flow for different hadron species in Au+Au collisions, color 
lines in red, blue and green show our predictions, for antiprotons, pions and 
kaons respectively, dashed lines experimental data from RHIC |48j . 



done retaining only the first term of the expansion in powers of — R. For 
Pt < 1.5 GeV/c, this term is small, however for larger px we should include 
additional terms. 

In fig. 3 we compare our results on the azimuthal angle dependence of 
the nuclear modified factor Raa{Pti ^V 9 ) a ^ different centralities with the 
PHENIX experimental data [49] . The model reproduces the shape for all 
centralities although the results are below the experimental data for low 
angles. 

In fig. 4 we plot v2/e part being e = J^l^ * 8 ^ the participant 

eccentricity as a function of the number of participants Na together the 
PHOBOS [50j[51j experimental data in Cu-Cu and Au-Au at y/s = 200 
GeV. It is observed good agreement. 

In fig. 5 we plot fa/e, being the participant eccentricity as a function of 
the ratio ^jj^jj where S is the surface of the collision for a given centrality. 
Together our results we show the experimental data at AGS, SPS and RHIC 
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Figure 3: Angular suppression for pions at different centrality, (10 — 20%, 20 — 
30%, 30 - 40%, 50 - 60%) results compared to PHENIX data. 
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Figure 4: Elliptic flow v2/e part as a function of the number of participants. 
Pink triangles and round points are the model results for Au-Au and Cu- Cu 
collisions at y/s = 200 GeV. 
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Figure 5: Scaling predictions for Au+Au at s/s = 200 GeV (stars). 



energies for different type of collisions [52]. A good agreement is observed. 
The scaling, experimentally observed, is also expected in our model. In fact, 
from equation (2) we obtain. 

which is proportional to J~p at high densities. On the other hand, from 
equations (22) and (20) we have 

v 2 (pt) = A(- [ 2 rf0 COS (20)(^ - 1)) (21) 

with 

e-' - F{-pf y\ 1 

F{p) 2^1 + W { ] 

As the bracket of the equation(24) is related to the centrality and A 

1 dN ch 
S dy 



depends only on p, i.e. depends only on i— f^. The equation (21) gives 



approximately the scaling law observed. 
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Figure 6: Dependence on the string density of the function described in the 
text 



In fig. 6 we plot the function — (e _p — F(p) 2 )/ F(p) as a function of 
p. We observe that increases up to p ~ 2.7 where it starts to decrease 
slowly. Therefore we predict this behavior for the elliptic flow. As the value 
p = 2.5 only is reached for central Au-Au collisions at RHIC or semicentral 
at LHC energies, only for this type of collisions will be sightly smaller at LHC 
energies. As the antibaryons(baryons) probe a higher string density than the 
mesons the elliptic flow for them will decrease for a less degree of centrality. 

In fig. 7 we plot V2/nq as a function of the kinetic energy over the con- 
stituent quarks of the observed particles. The recombination models [12] [J3] 
predicted an scaling law for intermediate transverse momentum, indepen- 
dently of mesons or baryons, all particles should lay in the same universal 
curve. We show our results for pions, kaons and antiprotons together with 
the experimental data [H]. It is seen that our model satisfies the scaling 
law as it was expected. In fact, in the percolation of strings, in addition 
to the mass effect due to a higher tension of the cluster which gives rise to 
enhancement of heavier particles, there is a flavouring effect related to the 
number of strings of the cluster, each one with flavour, which gives rise to 
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Figure 7: The red, green and blue solid lines stand for the model results for 
antiprotons, kaons and pions respectively compared with the experimental 
data (errors are smaller than symbols used). 



a higher probability to combine them, producing more baryons |39j [40] . In 
order to obtain close formula we incorporate this, only in approximate way 
as was explained before, however our results for pions, kaons and antipro- 
tons are laying approximately in the same curve as the experimental data, 
for KeT/ng < 0.2 GeV, corresponding to low pt- For higher p? the result is 
below the experimental data. 

In figs. 8, 9, 10 we show our results at yfs = 5.5 TeV compared with s/s = 
200 GeV for different hadrons. For the three particles we obtain a moderate 
increase of elliptic flow. This is a consequence that the string density at 
both energies for minimum bias Au-Au collisions is below 2.7. Above this 
density, the elliptic flow decreases what happened for Pb-Pb central collisions 
at LHC. 

In fig. 11 we show our prediction for t> 2 for pp collisions at ^/s = 14TeV 
for pions and kaons. A sizable elliptic flow is obtained, this result is not 



17 




Figure 8: Expected elliptic flow for antiprotons by our model, green lines 
RHIC and red lines LHC energies. 
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Figure 9: Expected elliptic flow pions by our model, green lines RHIC and 
red lines LHC energies. 
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surprising, given the high density partonic structure obtained in pp at LHC 
which is similar, see table 1, to Cu-Cu collisions The high density partonic 
structure of pp at LHC can give rise to other effects observed in AA collisions 
at RHIC as the high particles and the existence of large rapidity long range 
correlation even at very large intervals of rapidity [53j[54|[55j. 

5 Conclusions. 

In the framework of the clustering of color sources, the elliptic flow, t>2, and 
the dependence of the nuclear modified factor on the azimuthal angle have 
been evaluated. 

The comparison with RHIC experimental data for < 1.5 GeV/c shows 
a reasonable agreement. In particular, the observed mass ordering is repro- 
duced. In the model, is obtained an analytical expression for the observed 
scaling of the elliptic flow normalized to the participant eccentricity on the ra- 
tio The analytical expression indicates that above a determined density 
p = 2.7 elliptical flow decreases slowly. This density is reached in Au-Au cen- 
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tral collisions both at RHIC and LHC energies, therefore we predict slightly 
decrease in this case. For Au-Au minimum bias the density is lower than 2.7 
and we predict a slightly increase in that case. 

We have computed the dependence of the nuclear modify factor on the 
azimuthal angle comparing our results with experimental data for different 
centralities. We reproduce the general trend of the data, although we are 
15% below for small angles. 

We predict sizable elliptic flow in proton proton collisions, consequence 
of the high partonic density reached at LHC energy. The main ingredient of 
the model is the interaction among strings or equivalently among the partons 
of projectile and target located at the end of each string. In this way, we 
have a high density partonic interacting medium. The reasonable agreement 
obtained in our model with the experimental data for < 1-5 GcV/c, 
confirms the common belief that elliptic flow reveals a strongly interacting 
medium created in the first stage of the collision. Our predictions will be 
tested with the incoming LHC experiments. 





N s 


P 


Hp) 


dN 

dy 


pp (m.b) 


8 


0.4 


3.6 


2.5 




14 


0.7 


3.4 


5.5 


Cu-Cu (m.b) 


100 


0.4 


3.6 


50 




170 


0.7 


3.4 


120 


Au-Au (m.b) 


390 


1.0 


3.6 


180 




550 


1.6 


3.6 


400 


Au- Au (0 - 10)% 


1600 


2.7 


4.0 


650 




2600 


4.8 


4.1 


1500 



Table 1: Values of iV s , p, k(p) and for pp, Cu-Cu, Au-Au minimum bias 

(m.b) and — 10% central Au-Au collisions. For each collision the above 
numbers correspond to yfs = 200 GeV and the below ones to yfs = 5.5 TeV. 
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